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1. Introduction 

Electrostatic forces - the strongest interparticle forces outside of the nucleus - account for 
the stability of matter over a broad range of length scales, from atoms to macromolecules to 
the myriad materials that surround us. Colloidal particles (nanometers to microns in size) 
and polymers can become charged in a polar solvent (e.g., water) through dissociation of 
counterions [ij. Repulsive Coulomb interactions between ions can then stabilize a suspension 
or solution against aggregation due to ubiquitous van der Waals attractive interactions [2]. 
Electrostatic interactions between ions largely govern the equilibrium thermodynamic and 
dynamical properties of charge-stabilized colloidal suspensions and polyelectrolyte solutions. 
Controlling mechanical and thermal stability is essential to many applications - from foods 
and pharmaceuticals to filters and photonic materials. 

Predicting the properties of such complex, multicomponent systems with accuracy 
sufficient to guide and interpret experiments requires realistic modeling of the interparticle 
interactions and collective behavior of many-particle systems. While the fundamental 
interactions are simple, the sheer number of particles and the broad ranges of length and time 
scales confront the modeler with significant computational challenges. A general strategy for 
mitigating such challenges is to "coarse grain" or "integrate out" the degrees of freedom of 
some components, reducing the original multicomponent model to a simpler model of fewer 
components |3l IH |5l [6] . The trade-off for so reducing complexity is that the simpler model 
is governed by modified (effective) interparticle interactions. 

This chapter is a "how-to" manual of sorts for implementing coarse-graining methods 
to derive effective interactions. For simplicity, we focus here on charge-stabilized colloidal 
suspensions. The same basic concepts apply, however, to a wide variety of soft (and 
hard) materials. To set the stage, we begin by defining the primitive and one-component 
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models and outlining the Poisson-Boltzmann (PB) theory of charged colloids. After 
reviewing the well-known cell model implementation of PB theory, we turn to an alternative 
implementation, based on perturbation theory, and derive microion distributions around 
colloids and effective electrostatic interactions. After some effort, a happy result emerges: 
the effective interactions remain relatively simple, at least in systems with monovalent 
microions. This fortunate outcome provides the foundation for further theoretical and 
simulation modeling to explore and facilitate design of novel materials. Finally, we peer 
over the horizon at the outlook for possible future research directions in the field. 

2. Primitive Model 

The primitive model of charged colloids and polyelectrolytes [H [2] idealizes the solvent as 
a homogeneous dielectric continuum of relative permittivity e. Dispersed throughout the 
solvent are macroions and microions, modeled here as charged hard spheres of radius a and 
valence Z (charge —Ze) and point ions of valence z. In a closed suspension, all particles are 
confined to the same volume V. In Donnan equilibrium, only the macroions are confined, 
while the microions can exchange (e.g., across a semi-permeable membrane) with a salt 
reservoir, here assumed to be a 1:1 electrolyte of monovalent ions. 



If sufficiently dilute, the reservoir can be reasonably modeled as an ideal gas. The 
reservoir fixes the chemical potential of salt in the suspension, /ig = 2kBT\n{noA^), where 
no is the number density of ion pairs, T the absolute temperature, and A the thermal 
wavelength, defining the arbitrary zero of the chemical potential. A bulk suspension of Nm 
macroions, Nc counterions, and A^^ dissociated pairs of oppositely charged salt ions contains 
A^_i_ = Nc + Ng positive and A^_ = Ng negative microions, for a total of A"^ = + 2Ns 
microions. Global electroneutrality constrains the macroion and counterion numbers by the 
condition ZN^ = zN^. 




Figure 1. Primitive model of charged colloids: spherical macroions (radius a, valence 
Z) and point microions dispersed in a dielectric continuum (left). Colloidal suspension in 
Donnan equilibrium across semi-permeable membrane with electrolyte reservoir (right). 



Coarse- Grained Modeling of Charged Colloids 3 

Accounting for excluded-volume and electrostatic (Coulomb) pairwise interparticle 
interactions, the Hamiltonian separates naturally into three terms: 

H = Hm + + Hmf_i ■ (1) 

The first term on the right side is the macroion Hamiltonian 

Hm = Hhs + ^ ^ Vmm{rij) , (2) 

where i^hs is the hard-sphere Hamiltonian (including kinetic energy), Vmmij'ij) = Z'^Xs/rij is 
the Coulomb pair potential between a pair of macro ions with center-center separation rij, and 
Xb = e'^/leksT) defines the Bjerrum length, the distance at which the Coulomb interaction 
energy between a pair of monovalent ions rivals the thermal energy. To simplify notation, 
in this chapter, Hamiltonians, pairpotentials, and all other quantities having dimensions of 
energy are expressed in thermal (fc^T) units. The remaining two terms on the right side of 
Eq. ([T]) are the microion Hamiltonian 

h.^k.+^j:^, (3) 

with kinetic energy and microion valences Zi = ±1, and the macroion-microion interaction 
energy 

Hm,=zxsj:j:^. (4) 



i=i j=i 



3. One-Component Model: Effective Hamiltonian 

Significant concentrations of salt pose severe computational challenges for large-scale 
simulations of bulk suspensions [3 [HI |9]. Consider that a suspension of macroions of radius 
a = 10 nm at 1% volume fraction and just 1 mM salt concentration contains, per macroion, 
O{10^) particles, all interacting via long-range Coulomb forces. Therefore, salt-dominated 
suspensions usually are modeled by first mapping the mixture onto a one-component model 
(OCM). 

For a suspension in Donnan equilibrium, this mapping (or coarse graining) operates on 
the semigrand partition function: 

Z = {{expi-H)),)m , (5) 

where ( )^ denotes a grand canonical trace over microion coordinates and ( )m a canonical 
trace over macroion coordinates. The partition function may be formally expressed in the 
form 



Z = (exp(-iJeff)) 



(6) 
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where 

HeS = Hm + ^fj, (7) 

is an effective Hamiltonian of a one-component system of pseudo-macroions and 

^f, = - ln(exp(-iJ^ - Hmi,))i, (8) 

is tlie grand potential functional of the microions in the presence of fixed macroions. Practical 
applications of the OCM require approximating Q^. For this purpose, Poisson-Boltzmann 
(PB) theory [H [21 [TOl [11] is a powerful approach. Below we briefly review PB theory and 
two common implementations: the cell model and the effective-interaction model. 



4. Poisson-Boltzmann Theory 

Poisson-Boltzmann theory is most elegantly formulated within the framework of classical 
density-functional theory of nonuniform fluids [121 [131 El ESI IIS]- Corresponding to the 
primitive model Hamiltonian [Eqs. ([2|)-(|1])], there exists a Helmholtz free energy functional 
F[nm{r),n±{r)], which (for a given external potential) is a unique functional of the macroion 
and microion number density profiles, ^^(r) and n-|-(r), varying with position r [131- This 
free energy functional separates, according to F = Fid + -^ex + Fext, into a (purely entropic) 
ideal-gas free energy functional of all ions. 

Fid = / dr 5^ n,(r){ln[n,(r)A=^] - 1} , (9) 

i=m,± 

an excess free energy functional, Fex = Fhs-I-Fei, due to hard-sphere (hs) and electrostatic (el) 
interparticle interactions, and a contribution Fgxt due to an external potential. Neglecting 
interparticle correlations (mean-field approximation), the electrostatic part of the excess free 
energy functional may be approximated as 

Fei = lJ drp(r)^(r) , (10) 

where p(r) = e[n+(r) — n_{r) — nf{r)] is the total charge density, including the negative 
charge fixed on the macroion surfaces of number density ^/(r), and 

^(r) = ^fdr'^^ (11) 



e J |r — r I 

is the reduced electrostatic potential at position r due to all ions. The reduced potential and 
total ion density are related via the Poisson equation, which may be expressed in the form 

VV(r) = -47rABK(r) - n_(r)] ; V^|surface = Z\B/a^ , (12) 

where the macroion charges are absorbed into a boundary condition at the macroion surfaces 
and the microion densities implicitly vanish inside the macroion cores. 

In a given external potential, the equilibrium densities of all ions minimize the total 
grand potential functional of the system p!3]. Alternatively, fixing the macroions and 
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regarding their charges as the source of the external potential, the equilibrium microion 
densities alone minimize the microion grand potential functional 

^^,[n±{r)] = F^[n±(r)] - j drn+(r) - 1^- j drn_(r) , (13) 

a Legendre transform of the microion free energy functional 



F,[n±{r)] = F^,id[n±(r)] + 1 ^ drp(r)xl/(r) 



(14) 



where 



i^^,id[n±(r)] = fdr ^ n,(r){lnK(r)A=^] - 1} (15) 

i=± 

is the ideal-gas free energy functional and the microion (electro) chemical potentials n± 
are identified as the Legendre variables. Note that depends parametrically on the 
macroion coordinates and that includes macroion-macroion Coulomb interactions for 
electroneutrality. Under the assumption that either the electrostatic potential or the electric 
field vanishes everywhere on the boundary of the volume V, the microion free energy 
functional also may be expressed in the form 

Ff.[n±ir)] = jdrJ2 ni(r){lnK(r)A=^] - 1} + j dr\V^\\ (16) 

Minimizing f2^[n±(r)] with respect to n±{r) now yields the Boltzmann approximation for 
the equilibrium microion densities 

n±(r) = n^"* exp [=!='?/' (r)] (fixed macroions) , (17) 

where the reference densities, n^'' = A~^exp(/i±), are the microion densities at the reference 
potential = 0. The microion grand potential is the value of the grand potential functional 
[Eq. ([T5]) ] evaluated at the equilibrium density profiles [Eq. f lT7|) ]: 



— J dv [n+(r) + n^{r)] — - J dr [r7,_|_(r) — r7,_(r) + r2/(r)]^/'(r) 



For a closed suspension (fixed particle numbers), the chemical potentials of the two 
microion species differ because of asymmetric interactions with the macroions: /i^ ^ fi-. 
Correspondingly, the reference densities also differ: ^ In Donnan equilibrium, 

however, exchange of microions with a salt reservoir shifts the intrinsic microion chemical 
potentials, /i™ = [5F^/5n±(r)]cq, by the Donnan potential ipD- 

/i± = fj}^ ±iJD = ln(noA=^) . (19) 

The total chemical potentials, and so too the reference densities, of the two microion species 
are thus equalized. The equilibrium microion density profiles are then given by 

n±(r) = noexp[=F^(r)] . (20) 
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The Donnan potential is interpreted physically as the change in electrostatic potential across 
the reservoir-suspension interface, and mathematically as a Lagrange multiplier for the 
constraint of global electroneutrality. 

Combining the Poisson equation for the potential [Eq. (fT2|) ] with the Boltzmann 
approximation for the microion densities [Eq. (JTTll or (1201) ]. the Poisson-Boltzmann equation 
becomes 

V^V^(r) = /t^sinh^(r) ; | surface = ^A^/a^ , (21) 

where kq = a/SttXb^ is the screening constant in the reservoir. Note that Eq. ( 12T|) is highly 
nonlinear. In the case of weak electrostatic potentials {ip <^ 1), the right side of Eq. ( 12T|) may 
be expanded in powers of ip. Truncating at linear order yields the linearized PB equation: 

V'V'(r) = ; VV^|surface = Z\B/a' . (22) 

Beyond the boundary condition at the macroion surfaces, the boundary-value problem is 
fully specified only by imposing another condition at the outer boundary of the system, 
which depends on the practical implementation of the theory. 

5. Cell-Model Implementation 




Figure 2. Cell model: a single macroion in a spherical cell. 

The anisotropic boundary conditions on the nonlinear PB equation [Eq. (!2T]) ] imposed 
by an arbitrary configuration of macroions render a general solution computationally 
daunting. In recent years, powerful ab initio methods have been developed for combining 
PB theory of microion density profiles with molecular dynamics [151 UHl [13 HE] or Brownian 
dynamics [191 ISO] simulation to evolve macroion coordinates according to derived forces. 
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Despite such advances, most applications of PB theory have been implemented within a cell 
model to facilitate numerical solution. 

In a seemingly bold reduction, the cell model represents a bulk suspension by a single 
macroion, neutralizing counterions, and salt ions confined to a cell of the same shape as 
the macroion (see Fig. |2]). Microion- induced correlations between macroions are simply 
ignored. For spherical colloids, the natural choice is a spherical cell centered on the macroion. 
Gauss's law then dictates that the electric field must vanish everywhere on the boundary 
of the electroneutral cell. With the potential and microion densities depending on only 
the radial coordinate r, the PB equation reduces to an ordinary differential equation for 
ilj{r) with boundary conditions ip'ia) = ZXb/o^ and ip'{R) = 0, where the cell radius R is 
commensurate with the average macroion density: = N^/V = 3/(47ri?^). For a closed 
suspension, the arbitrary location of the reference point of the electrostatic potential (where 
ip = 0) is usually chosen as the cell boundary: ip{R) = 0. In Donnan equilibrium, the 
potential is conventionally chosen to vanish in the reservoir, in which case the boundary 
value of the electrostatic potential is identified as the Donnan potential: ipi^R) = ipD 7^ 0. 

An appealing feature of the cell model is the simple analytic relation between the bulk 
pressure p (in thermal units) and the microion densities at the cell boundary: 



Although first derived within the mean- field PB framework [21], this pressure theorem proves 
to be exact within the cell model p^, i.e., valid also for correlated microions. In Donnan 
equilibrium with an ideal-gas reservoir, the osmotic pressure 11, i.e., the difference in pressure 
between suspension and reservoir, is then given by 



Within PB theory, the osmotic pressure is strictly positive [23], following from Eqs. fl20l) and 
( l24l) and the inequality, coshx > 1. 

The cell model provides one means of implementing PB theory by approximating the 
microion grand potential [Eq. ([71)] in the one-component mapping of the primitive model 
(Sec. [3]). Reducing a suspension to a single macroion in a cell with isotropic boundary 
conditions facilitates solution of the nonlinear PB equation. The cost of incorporating 
nonlinear microion screening, however, is neglect of correlations between macroions. 

6. Effective-Interaction Implementation 

An alternative implementation of PB theory, based on the one-component mapping, 
focuses on effective interactions derived from perturbative expansion of the microion grand 
potential about a reference system, namely, a uniform plasma of microions unperturbed 
by the macroions. By incorporating macroion interactions, this approach can model both 
thermodynamic and structural properties of colloidal suspensions. For reviews of effective 
interactions in colloidal suspensions, see refs. [31 HI El E]. Several statistical mechanical 



p = n+{R) +n_{R) . 



(23) 



U = n+{R)+n^{R)-2no . 



(24) 
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frameworks have been developed. Density-functional theories [211 E51 EZ] expand the 
ideal-gas part of the microion free energy functional [Eq. f lT5|l ] in a Taylor series in powers 
of deviations of the microion density profiles from their mean values n±: 



i^M,id[^±(r)] = {N,[Hn^^') - 1] + y"dr 



2ni 



+ ■■■] ■ (25) 



i=± 

Distributon function theories are based on extensions of the Debye-Hiickel theory of 
electrolytes j28j. Response theories [291 EOl Ell |32] are based on a similar perturbative 
expansion of the microion grand potential functional [Eq. (|8|)]: 

n^ = no+ [ dX {Hmf,)x , (26) 

where Qq = — ln(exp(— if^))^ is the grand potential of a reference suspension in the 
absence of an external colloidal potential and {Hm^)x denotes an ensemble average of the 
macroion-microion interaction energy in a system in which the macroions are "charged" to a 
fraction A of their full charge, which can be related to the macroion-microion pair potentials 
Vm±{r) = ^ZXb/t and the densities, ^^(r) and ^±(r), of macroions and microions: 

{Hm±)x = j dr j dr' Vm±{\Y -v'\)nm{Y) {n±{v'))-^ . (27) 

Expanding the microion densities about a reference microion plasma in powers of the 
macroion "external" potential, 

0±(r) = j dT'vm±{\v - r'\)nm{v') , (28) 

and truncating the series at linear order yields the linear-response approximation 

n±(r) = Y,I dr'x±.(|r-r'|)0,(r') . (29) 

i=± 

The linear-response functions 

%(k-r'|)=(##) (30) 



describe the response of the reference plasma to the macroions and are related to the plasma 
pair correlation functions hij{r) via [33] 

Xij{r) = -ni[SijS{r) + njhij^r)] . (31) 

The neglected higher-order terms in Eq. ( l29l) involve nonlinear response functions and many- 
particle correlations. 

Approximating the pair correlation functions of the uniform plasma by their asymptotic 
(long-range) limits yields the so-called random phase approximation 

g-/tr 

hij{r) = -ZiZjXs—^ , (32) 
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where k = y^ATrXsn^ is the screening constant in the system, which differs from that in the 
reservoir kq [of. Eq. fl2T]) ]. For consistency, here represents the total density of microions 
in the free volume, i.e., the volume not excluded by the macroion hard cores. 

Inserting the linearized microion densities [Eq. ( 129|) ] into Eq. ( 127|) recasts the effective 
Hamiltonian (microion grand potential) in the form of a sum of effective interactions: 

^ Nm 

Hes = + Hhs + - ^ Ves{rij) , (33) 

where the volume energy E^j is the microion grand potential for a single macroion, feff(r) is 
an effective electrostatic pair potential between macroions, and neglected higher-order terms 
involve sums over effective many-body interactions. The volume energy takes the general 
form 

E, = no + ^N^Vin^iO) + i(iV+ - N_)^n , (34) 

where the first term on the right side accounts for the microion entropy, the second term 
for the macroion-microion interaction energy, and the last term is the Donnan potential 
energy resulting from the electroneutrality constraint. The effective macro ion-macroion pair 
potential can be expressed as 

Ves{r) = V.rnrn{r) + t^ind('^) , (35) 

which comprises the bare Coulomb pair potential Vmmi'f') and a microion-induced pair 
potential 

Vind{r) = / dr' [n+(r') - n_(r')]f„+(|r - r'|) . (36) 



To evaluate find('"); it is helpful first to consider the isotropic microion density profiles around 
a single macroion, which are obtained from Eqs. (!29|) - (!32|) : 



".W^E/'>r'x-..(l-r'IK..(/) 



= ±n±(^nf,XB J dr' ^—j-Vm+{\r - r'\) - Vrn+{r)^ . (37) 

To ensure exclusion of microions from the macroion hard cores, the macroion-microion pair 
potentials can be extended inside the core: 

{1/r , r > a 
(38) 
a/a , r < a , 

where a is a constant to be chosen such that n±{r) = for r < a. After evaluating the 
convolution integral 

e'^'^ e""" 
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with 

g{r,r') = j dfi f{\r - r'\) ; fi = r ■ r' /{rr') , (40) 

one finds [from Eqs. ( 1371) and ( l38l) ] a = Ka/{1 + na) and 

— , r > a . (41) 

1 + /ta/ r 

Substituting Eq. f HT]) into Eqs. f l5^ and f l5B]) . and evaluating tlie integral (for r > 2a) 

" dr' r'e-''^'^7(r, r') = ^ f - ^^e'^^^ , (42) 

yields a screened- Coulomb (Yukawa) effective pair potential: 

, r > 2a . (43) 

1 + /ta / r 

Finally, from Eqs. (|36l) and (|4T1) . we liave 

^ind(O) = , (44) 

a + K ^ 

whicfi allows the volume energy to be written more explicitly: 

£„.n„-iiV„^-l(^^^. (45) 

" 2 a + K-1 2 A/^^ ^ ' 



In passing, we recall that Eq. ( 1431) is the basis of the classic Derjaguin-Landau-Verwey- 
Overbeek (DLVO) theory [3H [35] of charged colloids, where it was first derived within a 
Debye-Hiickel approximation, without the associated volume energy and without excluded- 
volume corrections. 

It is important to remember that the one-body volume energy, although independent 
of macroion coordinates, does depend on the average macroion density, and therefore 
can influence thermodynamic properties. Similarly, the effective pair potential is density- 
dependent. Implications of density-dependent effective interactions for thermodynamic 
stability and consistency have been discussed extensively in recent years [361 EH EHJ [391 

MM- 

Bulk thermodynamic and structural properties can be calculated ultimately by inputting 
the effective interactions, and Ves{r), into statistical mechanical theories or simulations 
of the OCM. Figures [3] and [H show sample results for the osmotic pressure of highly 
charged colloidal suspensions, calculated using methodologies described in refs. [HI [121 [13]. 
Predictions of the PB cell and effective interaction (one-component) models agree closely 
with simulations of the primitive model [7] (Fig. [3]) up to moderate electrostatic coupling 
strengths (P = A^/a < 1) and with experimental data [lH [15] (Fig. [1]) over a considerable 
range of colloid volume fractions. To achieve such agreement, the linear- response theory has 
been merged with charge renormalization theory to incorporate the important concept of an 
effective macroion valence [16]. By subsuming within the effective valence those counterions 



Coarse- Grained Modeling of Charged Colloids 11 



Volume Fraction (|)=0.01 




0.01 0.1 1 

Coupling Strength r='k^/a 

Figure 3. Pressure vs. electrostatic coupling constant F = As/a for salt-free suspensions 
with colloid volume fraction (j> = 0.01 and several bare valences Zg. Predictions of the 
effective interaction model (curves) are compared with primitive model simulation data 
(symbols) [7]. 

that are strongly associated with the macroions, the renormahzed theory includes most 
of the nonlinear response inherent in the PB equation. Comparably accurate results are 
obtained with a renormalized jellium theory [l?! HSl HH]- The effective interaction model also 
predicts radial distribution functions in close agreement with simulations of the primitive 
model HHIig. 

7. Outlook 

To illustrate the essence of coarse-grained modeling in a friendly context, the discussion in 
this chapter is limited to the primitive model of charged colloids, focusing on monodisperse 
suspensions of microspheres. For this system, the effective electrostatic interactions 
derived using the coarse-graining scheme described in Sec. |6] are relatively simple and lead 
to predictions for thermodynamic and structural properties that agree closely with the 
Poisson-Boltzmann cell model, detailed simulations of the primitive model, and experiment. 
Tremendous opportunities now lie ahead for extending the general methods outlined here and 
applying them to more complex systems, such as mixtures of colloids differing in size, shape, 
and charge, to suspensions of particles with anisotropic (patchy) charge distributions (e.g.. 
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a=10nm, Z=628, c=10mM 




Volume Fraction log^Q(|) 

Figure 4. Osmotic pressure vs. colloid volume fraction for suspensions with macroion 
radius a = 10 nm, bare valence Zq — 628, and salt concentration Cg = 10 niM. Predictions 
of the effective interaction model (solid curve) and the PB cell model (dashed curve) are 
compared with experimental data (symbols) [44, 45]. 

Janus particles), as well as to different interparticle interactions, such as dipolar interactions. 
As computing power grows, the demand for coarse-grained models will likely persist, to guide 
and interpret simulations and to facilitate exploration of increasingly complex materials. 
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